Domain Coarsening in 2-d Ising Model: Finite-Size Scaling for Conserved Dynamics 
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We quantify the effect of system size in the kinetics of domain growth in Ising model with 50:50 
composition in two spatial dimensions. Our estimate of the exponent, a = 0.334 ± 0.004, for the 
power law growth of linear domain size, from Monte Carlo simulation using small systems of linear 
dimensions L = 16, 32, 64, and 128, is in excellent agreement with the prediction of Lifshitz-Slyozov 
(LS) theory, a — 1/3. We find that the LS exponent sets in very early and continues to be true 
until average size of domains reaches three quarters of equilibrium limit. 



A homogeneous binary mixture, A+B, becomes un- 
stable when quenched below the co-existence curve and 
moves towards the new equilibrium state via forma- 
tion and growth of domains rich either in A- or B- 
particles This domain coarsening is a scaling phe- 
nomena, e.g., the two-point equal time correlation func- 
tion, C(r,t), which characterizes the domain morphol- 
ogy and growth, exhibits the scaling behavior Q 



finite £(t) as 



C(r,t)=C(r/£(t)). 



(1) 



In Eq.([T]), C{x) is a scaling function independent of the 
average domain length £(t) which grows with time as 
£(t) ~ t a . Associating the rate of domain growth with 
the gradient of chemical potential, one can write [l[ 



d£(t) 
dt 



(2) 



where a is the A-B interfacial tension. Solving Eq.Q 
one gets a = 1/3, known as the Lifshitz-Slyozov (LS) 
3] growth law. 
While recent focus has been in more realistic systems 
understanding of kinetics of phase separation even 
in simple spin-i Ising model, where an up spin cor- 
responds to A-particle and a down spin to B-particle, 
appears to be incomplete. For conserved order param- 
eter, where composition of A- and B-species remains 
fixed during the entire evolution, dynamics in the Ising 
model is implemented via Kawasaki exchange mecha- 
nism 0, Q where for a Monte Carlo (MC) move inter- 
change of positions between a pair of nearest neighbor 
(nn) spins is tried. Earlier studies Q of phase ordering 
in conserved Ising model with 50:50 composition, most 
of which were based on MC simulation for very short 
period, reported estimates of a £ [0.17,0.25], deviating 
drastically from expected LS law. Even arguments in 
favor of logarithmic growth were proposed [8j . 

To resolve the controversy, Huse @ argued for the 
need to add an additional term oc l/£(t) 3 , to Eq.©, ac- 
counting for an enhanced interface conductivity, which 
brings in a correction to the instantaneous exponent at 



£{t) 



(3) 



Indeed, consistency of growth exponent with the LS law 
was established for 50:50 binary mixture 10] as well as 



multi-component mixture 11], in the limit i(t —> oo) 
— > oo. Present work, however, shows that observation 
of LS value of the exponent only for large £{t) is due to 
the mixing of a time dependent bare length in £(t). 

Most of the studies, till date, stressed on using large 
systems, with the anticipation of a strong finite-size ef- 
fect combined with the expectation that the LS law will 
be realized only in the large £{t) limit, and, of course, to 
obtain better self-averaging. The objective of this paper 
is to study finite-size effect on domain coarsening in 2-d 
conserved Ising model and understand the behavior of 
en from appropriate scaling analysis [T3 j . First effort 
in this direction was by Heermann, Yixue and Binder 
fl3l j for extreme asymmetrical composition, however, 
was never followed up systematically. In this paper, 
we show via application of finite-size scaling method 
that the LS value of a sets in very early and effect of 
size is rather small so that using a system of size as 
small as L 2 = 32 2 , one can confirm the LS growth law 
unambiguously. Such knowledge will be useful for com- 
puter simulation of more realistic systems by avoiding 
unnecessary large systems, thereby accessing large time 
scale. We also observe that systems do not provide self- 
averaging proportionate to their sizes [3]. 

In Fig HI we present snapshots during the evolution of 
an Ising system, starting from a 50:50 random mixture 
of up and down spins, on a square lattice of linear size 
L = 128, obtained via MC simulation at temperature 
T = 0.6T C , T c being the critical temperature. The times 
at which the shots were taken are mentioned on the fig- 
ure in units of Monte Carlo steps (MCS) where each 
MCS consists of exchange trial over L 2 pair of spins. 
Periodic boundary conditions were applied in both x- 
and y- directions. While the last snapshot corresponds 
to a situation when A and B phases are completely sep- 
arated, the one at t = 4.5 x 10 6 MCS represents the 



2 



situation when finite-size effect began to enter, which 
will be clear from subsequent discussion. Note that all 
physical quantities were calculated from pure domain 
morphology after eliminating the thermal noise via a 
majority spin rule where a spin at a lattice site j was 
replaced by the sign of majority of the spins sitting at 
j and nn of j . 
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FIG. 1: Evolution snapshots of domain coarsening in 2-d 
conserved Ising model following a quench from high tem- 
perature random state to T = 0.6T C , for a system of size 
L 2 = 128 2 . The black dots mark the location of A-particles 
while the B-particles are unmarked. 

FigH shows the plot of £(t) vs t for L = 32,64, and 
128, where £{t) was calculated from the first moment of 
the domain- length distribution function P(£k) with 
length £k being the separation between two successive 
interfaces in x- or y- directions. For L = 32, the data 
were averaged over 2000 independent initial configura- 
tions, for L = 64 and 128, averaging were done for 1000 
and 40 initial configurations, respectively. The flat re- 
gions of the data sets correspond to the situation when 
the systems reached their final equilibrium states, thus 
domains cannot grow beyond this. This limiting value, 
for present method of calculation, comes out to be ~ 
L/2 = £ max . The last snapshot in Fig Q] represents such 
a situation. In the inset of FigJ5J we demonstrate the 
scaling behavior of the correlation function, as embod- 
ied in Eq.((T]), for L = 128. The data collapse upon di- 
viding r by £(t) is good starting from as early as t = 10 3 
MCS till t = 4.5 x 10 6 MCS when the finite-size effect 
begins. Apparently, as is clear from the plot of £(t) vs 
t as well as the scaling behavior of C (r, t) for very ex- 
tended time and length scale, the size effect is negligible 



almost upto f max . However, to make quantitative state- 
ment about the extent of finite-size effect, appropriate 
scaling analysis is called for. 
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FIG. 2: Plot of average domain size as a function of time for 
three different system sizes L = 32, 64, and 128, as indicated. 
The inset shows a scaling plot of C(r,t) vs r/£(t) for L — 
128, for three different times, as mentioned. 

In analogy with critical phenomena, one can con- 



struct the finite-size scaling ansatz 13|, 1 1 51 ] by identi- 
fying £(t) with the equilibrium correlation length £ and 
1/t with temperature deviation from the critical point. 
At this stage we would like to quantify the growth of 
£(t) as 



£{t) = £(t ) + At° 



(4) 



where, to £(to), we assign the meaning of initial domain 
size when the system becomes unstable at time to af- 
ter the quench and the measurement of t starts from to ■ 
Note that slightly poor data collapse in Fig[2]for the ear- 
liest time is primarily because C (r, t) should have been 
plotted as a function of r/[£(t) — £{to)], not r/£(t). But 
C(r,t) contains information about £{t) and subtracting 
the influence of £(to) is a rather challenging task. 

Eq.Q is valid in absence of any finite-size effect. 
However, when £(t) is comparable to £ max , finite-size 
effect comes in and Eq.(U]) needs to be modified by ac- 
counting for the size effect as 



£{t) - £(t ) = y(x)t° 



(5) 



In Eq.([5]), y{x) is a scaling function, independent of the 
system size, which depends upon the scaling variable 
x = (£ max - £(to))/t a . By closely observing Eqs. flD 
and ([5]), one can write down the following behavior of 



3 



y(x) in the limiting situations: 



lim 

X— >0(i— >oo;£ n 

lim 



k<00) 



y(x) = x, 



>oo) 



y{x) = A. 



(6) 



(7) 



In FigEl we plot y(x) = [£(t) - £(t )] / >t a as a function 
of x/(x + xq) with xq — 5, for which we have varied 
a and ^(fo) to get optimum collapse of data from dif- 
ferent system sizes which is obtained for the choices 
£(t = 20) ~ 3.6 and a = 0.334. Note that £(t ) in 
our analysis is a bare length independent of time and 
the scaling behavior (5) will be obtained when this is 
chosen appropriately. These numbers, as expected, pro- 
vide a constant value of y(x) in the region unaffected 
by finite system size, which should be identified with 
the growth amplitude A. The arrow in Fig|3] marks the 
location where y{x) starts deviating from its constant 
value, signaling the appearance of finite-size effect at 
£{t) = (0.75 ± 0.05)^ max - This value is substantially 
large compared to expectation. Note that the snapshot 
at 4.5 x 10 6 MCS in FigfT] corresponds to this length. 
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— Eq.(8) 
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FIG. 3: Finite-size scaling plot of y(x), with £(to = 20) ~ 
3.6 and a — 0.334, as a function of x/(x + xo); Xo = 5. 
The continuous curve is a fit to Eq.((8]) with the best fit 
parameters mentioned in the text. The arrow roughly marks 
the appearance of finite-size effect. Inset: \£(t) — £(to)]~ 3 vs 
t~ , for L = 64, where the continuous line is a straight one 
with slope 39. 

Considering the limiting behaviors ((6|) and ([7]). we 
construct the following functional form of y(x) as 



y(x) 



Ax 



x + l/(p + qxP)' 



(8) 



The continuous line in Figgis a fit to the form ([5]) with 
A ~ 0.294, p = 3, q = 6400, ft = 4, which has the con- 
vergence lmxr^oo y(x) w — fx~ n }; n — 5. Of course, 
an exponential behavior can also not be ruled out. This 
behavior may be compared with much slower conver- 
gence of such scaling function in dynamic critical phe- 
nomena 0. In the inset of FigH we plot [e(t)-e(t )}- 3 
vs i _1 for L — 64 on a log-scale to bring visibility to 
a wide range of data. The continuous line there is a 
plot of the form ax with a ~ 39 = 1 /A 3 . The linear 
behavior of data starting from very early time justifies 
the introduction of £(to) again. 
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FIG. 4: Plot of instantaneous exponent Oi as a function of 
l/£'{t) for three different choices of £ s , for L = 64. The 
dashed straight lines have slopes —1.19,0, and 0.49, respec- 
tively. The arrow on the ordinate marks the value a = 1/3. 
Inset: a, vs l/£'(t) for £ s = 3.6 and L = 16, 32, and 64. 

Next we introduce a length l s to write 

£'{i) = £{t) - 4 = [£(to) - £ s ] + At a , (9) 

and define an instantaneous exponent a, = 
d[]n£'(t)]/d[hLt] to get 



£(to) 



t'(t) 



(10) 



According to Eq.(|T0|). when a; is plotted as a func- 
tion of l/£'(t), asymptotically one expects linear be- 
havior with a y-intercept equal to a. FigU shows 
such plots for l s = 0.0, 3.6, and 5.0, as indicated. 
The dashed lines have y-intercept a — 1/3 and slopes 
to = — [£(to) — £ s ]/3. The consistency of actual data 
with the dashed lines is remarkable, particularly, the 
behavior of <X; for £ s = £(t = 20) ~ 3.6 again speaks 



4 



for the choice of £(to) and strongly indicates that the 
LS scaling behavior sets in very early. Also see the in- 
set where we show on vs l/£'(t) for £ s ~ 3.6 for varying 
system sizes L = 16, 32, and 64. To this end, it is worth 
mentioning, as suggested by Huse, a term oc l/£(t) 3 in 
Eq.([2]), could also be motivated by introducing curva- 
ture dependence in a as a/[l + 2S/£(t)], S being the 
Tolman length [17j |. For a negative value of 5, one then 
would obtain time-dependent a;, as embodied in Eq.Q. 
However, our observation of negligible correction to the 
exponent, starting from very early time, is consistent 
with the fact that Tolman length is absent in a sym- 
metrical model fljj . 

The appearance of growing oscillation in ati around 
the mean value was also pointed out by Shinozaki and 
Oono fl9j | . In a finite system, as time increases, for an 
extended period of time two large neighboring domains 
of same signs may not merge, thus lowering the value 
of at, after which when they meet brings in drastic en- 
hancement in on . This character is in fact visible in the 
direct plot of £(t) for L = 128 in Fig[2] This oscilla- 
tion could be a route to an error if one obtains a from 
least square fitting. However, apart from averaging over 
huge number of initial configurations, it could also be 
reduced by considering times well separated from each 
other while taking instantaneous derivative. 

In conclusion, this paper contains a comprehensive 
finite-size scaling analysis of domain coarsening in a 
phase separating system. Our accurate and appropri- 
ate estimate, for which we quote 0.334 ± 0.004, of the 
growth exponent is almost coincident with the expected 
LS value a — 1/3, within tiny error bar. As opposed to 
the earlier understanding, correction appears to be very 
small, thus LS scaling regime sets in very early. Very 
small primary finite-size effect is a welcome message 
which is suggestive of avoiding large systems, rather 
focusing on accessing long time scale. Even though 
size effect may be situation and system dependent, re- 
cent study [20| | of phase separation in a binary fluid 
provides qualitative agreement with the finding of the 
present work. Nevertheless, one should be prepared to 
encounter stronger size effect in more complicated situ- 
ations, e.g., systems giving anisotropic patterns 

In a future paper, we will address the similar is- 
sue for phase separation, in different spatial dimensions 
with varying compositions where larger curvature in do- 
main morphology may delay the convergence to the LS 
regime. It will be interesting to look at the temperature 
dependence, in particular, the behavior of the scaling 
function y(x) as one approaches the critical tempera- 
ture T c and thus the validity of scaling of C(r,t) for 
extended period. A comparison of the finite-size effect 
in non-conserved order parameter situation, e.g., order- 
ing in a ferromagnet, with the present one is also under 



investigation. 
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